%Compare P0requencies oP0 given species
species;
A = load('ssaminout.txt');
P0f = [];
P0f = [P0f; mean(A(:,P0))];
P0f = [P0f; mean(A(:,P0_TALA_KRAB))];
P0f = [P0f; mean(A(:,P0_TALA_VP16))];
P0f = [P0f; mean(A(:,P0_TALA_KRAB_KRAB))];
P0f = [P0f; mean(A(:,P0_TALA_VP16_KRAB))];
P0f = [P0f; mean(A(:,P0_TALA_VP16_VP16))];

P1f = [];
P1f = [P1f; mean(A(:,P1))];
P1f = [P1f; mean(A(:,P1_TALA_KRAB))];
P1f = [P1f; mean(A(:,P1_TALA_VP16))];
P1f = [P1f; mean(A(:,P1_TALA_KRAB_KRAB))];
P1f = [P1f; mean(A(:,P1_TALA_VP16_KRAB))];
P1f = [P1f; mean(A(:,P1_TALA_VP16_VP16))];

P2f = [];
P2f = [P2f; mean(A(:,P2))];
P2f = [P2f; mean(A(:,P2_TALB_KRAB))];
P2f = [P2f; mean(A(:,P2_TALB_VP16))];
P2f = [P2f; mean(A(:,P2_TALB_KRAB_KRAB))];
P2f = [P2f; mean(A(:,P2_TALB_VP16_KRAB))];
P2f = [P2f; mean(A(:,P2_TALB_VP16_VP16))];

P3f = [];
P3f = [P3f; mean(A(:,P3))];
P3f = [P3f; mean(A(:,P3_TALB_KRAB))];
P3f = [P3f; mean(A(:,P3_TALB_VP16))];
P3f = [P3f; mean(A(:,P3_TALB_KRAB_KRAB))];
P3f = [P3f; mean(A(:,P3_TALB_VP16_KRAB))];
P3f = [P3f; mean(A(:,P3_TALB_VP16_VP16))];

Pf = [P0f P1f P2f P3f];
Pf = Pf./max(max(Pf));